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Abstract 

The production of a hard and isolated photon accompanied by one or two jets 
in large-Q 2 deep inelastic ep scattering is calculated at next-to-leading order. We 
include consistently contributions from quark-to-photon fragmentation and study 
various differential cross sections and their dependence on isolation cut parameters. 
Numerical results relevant for HERA experiments are presented. 
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1 Introduction 



In the past, measurements of prompt photon production at both fixed-target facilities and 
hadron-hadron colliders, have extensively been used to constrain the gluon distribution of 
the proton JT|. Only recently the first data on prompt photon production in high-energy 
ep collisions have been reported f2|. Due to the presently still limited statistics the mea- 
surements are confined to prompt photons produced in photoproduction reactions, i.e. to 
ep collisions with almost real exchanged photons (Q 2 ~ 0). As is well-known, photopro- 
duction processes at high energies proceed by two distinct mechanisms. The incoming 
photon can couple either in a point-like manner to the incoming quark or antiquark (direct 
process) or hadronically as a source of quarks and gluons which in turn take part in the 
subsequent hard scattering process (resolved process). Therefore an important advantage 
of photoproduction measurements is to provide additional constraints on the quark and 
gluon content of the photon as suggested many years ago by Aurenche et al. [||]. 

By contrast, prompt photon production in deep inelastic scattering with large Q 2 , Q 2 <; 10 
GeV 2 , is fully determined by the direct process and does not need any non-perturbative 
input for the parton content of the photonQ as in photoproduction. Therefore this process 
is sensitive only to the parton distribution functions (PDF's) of the proton. The possible 
information on the proton PDF's would be complementary to the F% measurement from 
inclusive deep inelastic scattering, since up- and down-type quarks contribute with dif- 
ferent weights. Of course, the cross sections for ep —>■ e^X at large Q 2 are smaller than 
the corresponding cross sections for almost real photons; but with the larger luminosities 
planned at HERA rather accurate measurements of various differential cross sections for 
Q 2 > 10 GeV 2 seem feasible. 

In addition to the perturbative direct production, photons are also produced through the 
"fragmentation" of a hadronic jet into a single photon carrying a large fraction of the jet 
energy 0]. This long-distance process is described in terms of the quark-to-photon and 
gluon-to-photon fragmentation functions (FF's). 

In order to unambiguously identify the prompt photon signal from the hadronic back- 
ground it is necessary to apply isolation cuts in the experiment. This has the effect of 
reducing the cross section due to a suppression of the photon fragmentation contributions. 
On the other hand, it has the advantage of eliminating to a large extent the dependence 
on the photon fragmentation function, which is a non-perturbative input and must come 
from other experiments designed to measure them. 

To obtain reliable predictions it is necessary to calculate the ep — * e^X cross section in 
next-to-leading order (NLO) of the strong coupling constant a s as has been done for pp 
collisions ||, e + e~ annihilation || [7|, ||, as well as photoproduction ||. The corresponding 
NLO calculation of cross sections for prompt photon production in ep scattering with large 
Q 2 has not yet been done previously, neither for the technically simpler case of inclusive 

*At intermediate Q 2 < 10 GeV 2 there might be a significant resolved contribution which depends on 
the quark and gluon distribution of the virtual photon. 
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cross sections, i.e. without any photon isolation cut, nor for the case with isolation cuts. 
Applicability of perturbative QCD requires that the scattering process is characterized 
by a large transverse momentum, provided either by the momentum transfer Q 2 or a 
large transverse momentum of the hadronic final state. We consider only the case where 
both Q 2 is large and the hadronic final state is characterized by a large Pt- One specific 
possibility is to consider the case where in addition to the photon also one or more jets 
are observed in the final state. The detection of an additional jet may also help to identify 
the prompt photon events in the actual experiment. In leading order (LO) the photon is 
produced by the Compton process 7* + q — > 7 + q, where 7* is a photon of high virtuality 
emitted by the incoming electron. This partonic photon production process contributes 
to the 7 + (1 + l)-jet final state in ep scattering (the proton remnant jet being counted 
as jet as usual). In NLO the final states are 7 + (1 + 1)- and 7 + (2 + l)-jets. 

The first NLO calculations for this case were done by two of us and D. Michelsen |T0|| . 
This calculation was restricted to the case of not too large Q 2 where it is possible to 
neglect the exchange of a Z boson. Moreover, in this previous work the fragmentation 
contribution was not taken into account. Therefore, photon-quark collinear singularities 
could not be absorbed into the fragmentation functions. Instead, these singularities had 
been removed by explicit parton-level cutoffs. As a consequence, the result depended 
strongly on these parton-photon cutoffs, in particular for subprocesses with an incoming 
gluon. These cutoffs are difficult to control experimentally, where hadrons combined into 



jets are observed and not the partons needed to define the cutoffs. In subsequent work [11 



we included the fragmentation contribution thereby avoiding the need to use parton-level 
cutoffs. The isolation criteria, which limit the hadronic energy in the jet containing the 
photon, are thus physical, i.e. correspond to selection criteria in the experimental analysis. 
In a later paper |12| we studied the sum of the 7 + (1 + l)-jet and 7 + (2 + l)-jet cross 
sections as a function of the momentum fraction carried by a photon inside a jet. We 
observed that this special cross section is sensitive to the fragmentation contribution, in 
particular to the quark-to-photon FF. 

In [|TTJ, only a few observables have been calculated, as for example distributions with 
respect to the transverse momentum, p?, and the rapidity, rj, of the photon or the most 
energetic jet for one particular choice of the photon isolation cut. In this paper we take up 
the topic of this earlier work. Besides several other observables which are of interest for 
analyzing upcoming experimental data from HERA we shall present results for px and r\ 
distributions already considered in JTI] for different isolation cuts. We also study possible 
scale dependencies to estimate the reliability of our predictions. As in |TT| we use the 
7*p center-of-mass system to define the kinematic variables. For most of the cross section 
calculations a particular cone algorithm is applied to define the parton jets and to isolate 
the photon signal. For a few cases we shall also make use of the hr cluster algorithm. 

The plan of the paper is as follows. In section 2 a brief outline of the theoretical back- 
ground to the cross section calculations as well as the technique of the calculation are 
given. The results for the various observables are presented and discussed in section 3. 
Section 4 contains a summary and some concluding remarks. 
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2 Subprocesses Through Next-to-Leading Order 



2.1 Leading-Order Subprocesses 

In leading order, the production of photons in deep inelastic electron scattering is de- 
scribed by the quark (antiquark) subprocess (see Fig. |l]) 

e(pi) + q(p 3 ) -> e(p 2 ) + q(p 4 ) + 7(^5) (1) 

where the particle momenta are given in parentheses. The momentum of the incoming 
quark is a fraction £ of the proton momentum P, p^ = £P. The proton remnant r has the 
momentum p r = (1 — £)P. It hadronizes into the remnant jet so that the process ([I]) gives 
rise to 7 + (1 + l)-jet final states. In the virtual photon-proton center-of-mass system the 
hard photon recoils against the hard jet back-to-back. In a leading-logarithmic calculation, 
the effects of higher-order processes show up only via the use of the scale- violating parton 
distributions of the proton. The PDF's are calculated using collinear kinematics so that 
the event structure is the same as given by the lowest-order subprocesses, which are of 
order 0(a 2 ) Q. 



lAAAAA/ 7 



AAA/W * 

7* q 
a) b) 

Figure 1: Feynman diagrams for 7* + q — > q + 7. 

To remove photon production by incoming photons with small virtuality and to restrict to 
the case where the scattered electron e(j>2) is observed, one applies cuts on the usual deep 
inelastic scattering variables x, y, Q 2 as measured from the momentum of the scattered 
lepton. In particular we restrict the calculation to values of Q 2 > 10 GeV 2 ; however, 
since very large Q 2 values are not relevant at HERA we can neglect contributions from Z 
boson exchange. In addition, to have photons 7(^5) of sufficiently large energy we require 
an explicit cut on their transverse momentum. Finally, a cut on the invariant mass of the 
hadronic final state is also applied. 

Both leptons and quarks emit photons. The subset of diagrams where the photon is emit- 
ted from the initial or final state lepton (leptonic radiation) is explicitly gauge invariant 
and can be considered separately. Similarly, the contribution from diagrams with a pho- 
ton emitted from quark lines is called quarkonic radiation. In addition, there are also 



2 Here and in the following we do not count the extra factor a from the eej* vertex. 
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contributions from the interference of these two mechanisms. The emission of photons 
from leptons is described by pure QED and can be predicted with high reliability. There- 
fore, the contributions from leptonic radiation will be suppressed by cuts on the photon 
emission angle [fLOR . In our numerical evaluation we include the remaining background 
from leptonic radiation as well as the interference contribution. 

At lowest order, each parton is identified with a jet and the photon is automatically 
isolated from the quark jet by requiring a non-zero transverse momentum of the photon 
or jet in the j*p center-of-mass frame. Therefore the photon fragmentation is absent in 
this order. 



2.2 Subprocesses to Next-to-Leading Order 
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Figure 2: Examples of Feynman diagrams for higher-order processes: a) 7* + q 
b) 7* + 9 -> q + q + 1- 



Q+9 + 1, 



At next-to-leading order, processes with an additional gluon, either emitted into the final 
state or as incoming parton, have to be taken into account: 



e(pi) + q(p 3 ) -> e(p 2 ) + q(p4) + t(Ps) + 9{Pe) 



(2) 



e(pi) + gips) -> e(pa) + g(p4) + 7(Ps) + ?(pe) 



(3) 



where the momenta of the particles are again given in parentheses. Examples of diagrams 
for 7*g — > 97(7 and 7*^ — > q'jq are shown in Fig. In addition, virtual corrections 
(one-loop diagrams at 0(a s )) to the LO processes (0) have to be included. The complete 
matrix elements for the processes and (|3]) are given in [13|. The processes (0) and 



(II) contribute both to the 7 + (1 + l)-jet cross section, as well as to the cross section for 
7 + (2 + l)-jets in the final state. In the latter case each parton in the final state builds a 
jet on its own, whereas for 7 + (1 + l)-jets a pair of final state partons is experimentally 
unresolved. The criteria for combining two partons into one jet will be introduced later. 
The contributions (||) and @ as well as the virtual corrections to (P are of order 0(a 2 a s ). 
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2.3 Fragmentation Contributions 



In addition to the direct production described in the last two subsections, photons can also 
be produced through the fragmentation of a hadronic jet into a single photon carrying 
a large fraction of the jet energy |4j]. This long-distance process is described in terms 
of quark-to-photon and gluon-to-photon fragmentation functions which absorb collinear 
singularities present in the NLO direct contributions of section 2.2. The corresponding 
fragmentation processes (see Fig. |3|) are 

e(pi) + q(p 3 ) -> e(p 2 ) + q(p 4 ) + g(p 6 ), (4) 

e(pi) + g(p 3 ) -> e(p 2 ) + q{pi) + q(p e ). (5) 

These processes are of order aa s whereas the photon FF is formally of order a, so that 
the LO fragmentation contribution is formally of order 0(a 2 a s ), i.e. of the same or- 
der as the NLO direct contribution. The fragmentation photons, sometimes also called 
bremsstrahlung photons, are emitted predominantly along the direction of motion of the 
parent quark or gluon. Because of the pointlike nature of the photon-quark interaction, 
it is possible to calculate the leading-logarithmic behaviour of the photon FF, including 
the corrections due to additional gluon emissions. The resulting FF's are in fact of order 
0(a/a s ) since they possess a logarithmic growth coming from the integration over the 
momenta of unobserved partons. Therefore in the leading-logarithmic approximation the 
fragmentation contribution is obtained from 0(a/a s ) FF's convoluted with corresponding 
0(aa s ) cross sections for the two-body subprocesses ([|) and @. The resulting contribu- 
tion is thus of order 0(a 2 ), i.e. of the same order as the LO non-fragmentation process 
(H|). For this reason it is sometimes argued that the fragmentation contribution should 
be combined with the LO direct process to provide the full LO physical cross section. 
Consequently, the calculation of the full cross section up to NLO would then also require 
the computation of the NLO corrections to the fragmentation contributions. On the other 
hand, it is well-known, that the fragmentation contribution in LO depends strongly on the 
factorization scale /ip which, however, is cancelled to a large extent by the /independence 
of the NLO contribution to the non-fragmentation part. For this reason and also since for 




Figure 3: Feynman diagrams giving rise to fragmentation contributions: a) '~f* + q — > q + g, 

h)Y + g^q + q- 
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an isolated photon the fragmentation contribution is small we shall take it into account 
only in LO in the same way as in our previous work [JIT, |l^ . 

The signature of the fragmentation contribution in LO is a photon balanced by a jet on 
the opposite side of the event and accompanied by nearly collinear hadrons on the same 
side of the event. This means that this contribution has a similar event structure as the 
LO direct contribution. 



2.4 Calculational Details 

The calculation of the NLO corrections was performed with the help of the phase space 
slicing method using a slicing parameter defined in terms of invariant masses. With this 
method it is straightforward to introduce the photon isolation requirement as well as to 
implement a jet definition which separates 7 + (1 + l)-jet from 7 + (2 + l)-jet final states. 
Phase space slicing based on invariant masses is also used to separate the photon-quark 
collinearly singular regions, however, using another independent cutoff parameter. The 
technical steps to apply the phase space slicing method in the present case are described 
in the following. 

In the calculation of the contribution to the 7 + (1 + l)-jet cross section we encounter the 
well-known infrared singularities. They appear in those phase space regions where two 
partons are degenerate to one parton, i.e. when the gluon becomes soft or two partons 
become collinear to each other. The singularities are assigned either to the initial state 
(ISR) or to the final state (FSR). Contributions involving the product of an ISR and a 
FSR factor are separated by partial fractioning. The FSR singularities cancel against 
singularities from virtual corrections to the LO process ([I]) . For the ISR singularities, the 
cancellation is incomplete and the remaining singular contributions have to be factorized 
and absorbed into the renormalized PDF's of the proton. 

To carry out these steps, the singularities are isolated in an analytic calculation using 
dimensional regularization. Since the corresponding calculations are too difficult for the 
exact cross sections of the processes (0) and (||) an approximate solution is required. To 
achieve this, the phase space slicing JTJJ is used first to separate the singular regions in 
the 4-particle phase space. Then, in these regions the matrix elements are approximated 
by their most singular contributions. Only for these approximate expressions and only in 
the singular regions the calculation is performed analytically. The separation of singular 
regions is obtained by applying a slicing cut Dq to the scaled invariant masses yij, where 
Vij = (Pi +Pj) 2 /W /7 had and the mass of the hadronic final state H^ad is defined by W^ &A = 
(P + g— P5) 2 . The cut j/q must be chosen small enough so that terms of order 0(i/q) which 
are neglected in the singular approximation are so small that an accuracy of a few per cent 
is achievable for the final result. Outside the singular regions the integrations are done 
numerically without any approximation and with 4 space-time dimensions. Physical cross 
sections, as defined in the next section, are obtained by adding the contributions from 
singular and non-singular regions as well as the virtual contributions and subtracting the 
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remaining ISR collinear singularities. In the final results, the dependence on the slicing 
parameter y$ cancels. This means the cut-off yfi is purely technical. The independence on 
the slicing cut y$ has been checked by explicit calculation for some special photon plus jet 



cross sections in |T(J. Further details and the derivation of the two-body matrix elements 



in the singular region together with the cancellation of the soft and collinear poles can be 



found in 13, 



In addition, the squared matrix elements for the processes (fj) and (§) have photonic 
infrared and collinear singularities, i.e. singularities due to soft or collinear photons. Since 
we require the photon to be observed in the detector the infrared singularity can not occur. 
In the numerical analysis we will introduce this condition by requiring a minimum on the 
transverse momentum of the photon. This cut removes also all collinear singularities due 
to initial state radiation. 

Final state collinear singularities due to photons are present and are treated again with 
the help of the phase space slicing method in a similar way as the quark-gluon collinear 
contributions. The phase space slicing parameter used to treat the photonic singularities 
can be chosen independently and is denoted by y$. As before, it has to be chosen very 
small so that the matrix element can be approximated by its singular part. For the 
subprocess (0), the phase space slicing is described by the squared invariant masses j/45 = 
(j>4 + P5) 2 /W^ad • I n the gluon-initiated process (§) one has two singular regions which 
are controlled by the variables y 45 and y 56 , respectively. In the regions y 45 , y 56 > yl 
the cross section is evaluated numerically in the same way as in [|Kj where these cuts 
were introduced as physical isolation cuts on the photon with sufficiently large isolation 
parameter y$. In this work the cuts on y 45 and y$Q are only technical since we include the 
contribution to the matrix element also in the regions 7/45 < y$ and y<j§ < y^. In these 
regions the matrix elements are collinearly divergent. The singularities are regulated 
by dimensional regularization, allowing us to absorb their divergent parts into the bare 
photon FF to yield the renormalized photon FF denoted by D q ^ r For the process (|2]) 
this procedure results in a contribution of the following form: 

\M\ 2 rq ^ qg = \M\ 2 rq ^ qg ® D g ^(z). (6) 

The matrix element |M| 2 on the right-hand side of (§) is the matrix element for the 
process 7*g — > qg whose Feynman diagram is shown in Fig. [3]a. There exists a similar 
expression for the subprocess 7*g — > 'jqq (Fig. |3|b). The photon FF D q ^(z) in (|6|) is 



given by [15 



ZWz) = D q ^(z, /4) + ^ (P^(z) In Z{1 + zj . (7) 

D q ^ 1 (z, n 2 F ) in (0) stands for the non-perturbative FF describing the transition q — ► 7 
at the factorization scale This function will be specified in the next section. The 
second term in (0), if substituted in (0), is the finite part due to the collinear photon- 
quark (-antiquark) contribution to the matrix element \M\ 2 * q ^ iqg integrated in the region 
2/45 < Uo after absorption of the divergent part into the non-perturbative FF. 
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Again the parameter is only a parameter used in intermediate steps of the calculation, 
introduced to separate divergent from finite contributions; the ^-dependence m © i s 
canceled by the dependence of the numerically computed 7 + (1 + l)-jet cross section 
restricted to the region j/45 > y$. Since the corresponding contributions to the matrix 
element in (^) are calculated in the collinear approximation, the result is valid only up to 
terms of order 0(uq). This requires to choose a very small value for yj- m Ref- HH it 



has been explicitly shown that the sum of all terms for the photon plus jet cross section 
becomes independent of y^ when y 7 is chosen small enough. 

In (|7|), P^> is the LO quark-to-photon splitting function 

P£Hz) = 1 + { \~ Z? (8) 

and e q is the electric charge of quark q. The variable z denotes the fraction of the quark 
momentum carried away by the photon. If the photon is emitted from the final state 
quark with 4-momentum p' 4 = p 4 +ps, then z can be related to the invariants 2/35 and 7/34: 

z = Vs5_ = 2/35 ^ 

1/34' 1/34 + 2/35 ' 

The fragmentation contribution is proportional to the cross section for 7*g — > qg which is 
0(aa s ) and is well-known. It must be convoluted with the function in ([FD to obtain the 
contribution to the cross section for 7*g — > q^g at 0(a 2 a s ). Equivalent formulas are used 
to calculate the fragmentation contributions to the channel (0) and in the case where the 
quarks in the initial and final state are replaced by an antiquark in (fj). 



3 Results 

3.1 Kinematical Selection Cuts and Other Input 

The results for the cross sections which will be presented in the following subsections 
are obtained for energies and kinematical cuts appropriate for HERA experiments. The 
energies of the incoming electron (positron) and proton are E e = 27.5 GeV and Ep = 820 
GeV, respectively. The cuts on the DIS variables are chosen as in our previous works 



11, 12 



Q 2 > 10 GeV 2 , W had > 10 GeV , 

(10) 

10~ 4 < x < 0.5, 0.05 < y < 0.95 
where Q 2 = —q 2 and q is the electron momentum transfer, q = pi — p 2 as usual. 
To reduce the background from leptonic radiation [ITH] we require 



90° < Q 1 < 173°, 6 ie > 10° (11) 
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where # 7 is the emission angle of the photon measured with respect to the momentum of 
the incoming electron in the HERA laboratory frame. The cut on # 7e (= #25), the angle 
between the momenta of the photon and the scattered electron, suppresses radiation from 
the final state electron. In the case of photon plus jet cross sections, the photon and the 
hadron jets J are required to have minimal transverse momenta in the j*p center-of-mass 
system (i.e. in the rest system of p\ — p<i + P where the remnant has p?,r — 0), 

Pt, 7 > 5 GeV, p T:J > 6 GeV . (12) 

Note that an event is rejected if we do not find at least one jet with p^v above the cut in 
(|l2l). If there are two partons with p? t j > 6 GeV, the event has a chance to be treated as 
a 7 + (2 + l)-jet event. If, after trying to recombine partons into jets (see below), only one 
jet has pt,j > 6 GeV, the event contributes to the 7 + (1 + l)-jet class. In the latter case, 
it may happen that there is an additional parton not combined into a jet with px.i < 6 
GeV (in a JADE-like jet algorithm applied in the HERA laboratory frame, such a low-pr 
parton would be recombined with the remnant jet in most cases). 

Different values of minimal p^s for the photon and the jet have to be chosen in order to 
avoid the otherwise present infrared sensitivity of the NLO predictions [[ll|]. This point 



will be studied in detail later. For inclusive photon cross sections, i.e. in the case where 
we do not perform a jet analysis of the hadronic final state, we replace the second of the 
conditions in (fl2] ) by a cut on the sum of transverse energies of all final state partons 

Et= E |Pr,<|>6GeV. (13) 

The PDF's of the proton are taken from [|1| (MRST). Recent updates of available PDF 
parametrizations |T7], |TB| (MRS99, CTEQ5) do not lead to markedly different results 
(roughly a 2 % (4 %) increase of the total cross section for MRS99 (CTEQ5) with respect 
to MRST; a version with enhanced d/u ratio of MRS99 does not lead to observable 
differences in the range of x and Q 2 considered here). a s is calculated from the two- 
loop formula with the same A- value (Ajjg(nf = 4) = 300 MeV) as used in the MRST 
parametrization of the proton PDF. The scale in a s and the factorization scale are set 
equal to each other and fixed to the largest Pt,j, except when we present results with 
other scale choices. For completeness, we also mention that the slicing cuts have been 
fixed at y£ = 10" 4 and = 10" 5 . 

The dependence of the 7+jet cross sections on the choice of the photon FF has been 
studied earlier [P2| . In the present work we choose the FF of Bourhis et al. (19| (BFG). 



This FF has been compared to the ALEPH e e — > 7 + 1-jet cross section |20[ and also 



to the inclusive photon cross section measured by OPAL [21 ] . Both data sets agreed well 



with predictions based on the BFG parametrization [^1|, In [12[ we studied the cross 
section differential with respect to the fraction of momentum z 7 carried by a photon inside 
a jet for several other photon FF's besides the BFG parametrization. 

Photon isolation is implemented with the help of the cone isolation method similar to 
the one used in the ZEUS experiment for photon production with almost real photons 
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0. This method restricts the hadronic energy allowed in a cone around the jet axis of 
the jet containing the photon. The same method is used also to define jets emerging in 
the event sample of 7 + 2-parton-level jets when two partons are combined. In the j*p 
center-of-mass frame, two partons i and j are combined into a jet J, when they obey the 
cone constraints R^j < R and Rj t j < R, where 



Ri,j = yJ(Vi-Vj) 2 + (4>i-4>j) 2 - (14) 

t]j (= — lntan(0j/2)) and 0j are the rapidity (polar angle 6j) and azimuthal angle of the 
recombined jet which are obtained from the formulae 

PT,J = PT,i + PT,j , 



( v.±. J. jfc • IJ- LJ -jJ 

VJ= PTJ+PTJ ' ( 15 ) 



ViPT,i + VjPTj 

Pt,% + Ptj 

^ _ 4>iPT,i + 4>jPT,j 
PT,i + PT,j 

For most of the results we choose R = 1. If not, we shall state the value of R explicitly. The 
azimuthal angle is defined with respect to the scattering plane given by the momentum 
of the beam and the momentum of the scattered electron. The rapidity is always defined 
positive in the direction of the proton remnant. The photon is treated like any other 
parton in the recombination process, i.e. when in (|14]) % or j is the photon, then J is 
called the photon-jet. To qualify a jet as a photon-jet, we restrict the hadronic energy in 
the jet by 

£7 = = 1 — Chad > 1 — Chad = z cut- (16) 

PT i7 + PT,had 

Pr,7 and pr,had denote the transverse momenta of the photon and the parton producing 
hadrons in this jet. Defining pr i7 _j et = +PT,had, (0) is equivalent to the requirement 

PT,had < Chad Pr, 7 -jet, (17) 

i.e. the ratio of the total px due to other particles (partons) than the photon is required 
not to exceed e° ad . For our predictions we shall choose different values for e° ad = (1 — z cut ). 
In this parameter was set equal to e{J ad = 0.11. 

It is known that the cone algorithm is ambiguous for final states with more than three 



particles or partons \22\. Since we have maximally three partons in the final state this 
problem is not relevant in our case. However, in some cases it may happen that two 
partons i and j qualify both as two individual jets % and j, or as a combined jet J. In 
these exceptional cases we count only the combined jet J to avoid double counting. 

The cone algorithm is problematic in experimental analyses due to its seed-finding mech- 
anism and due to overlapping cones. These problems are avoided with the kx jet finding 



algorithm [23 . In our theoretical calculations, the kx algorithm can be incorporated quite 
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easily: partons i and j (where one of them may be the photon) are combined if they fulfill 
the condition 

Rij < R with Rij = yJ( Vi - Vj y + (0, - 0,) 2 - (18) 

The resulting kinematic variables of the combined jet are calculated with the same for- 
mulae (p~5|) as in the cone algorithm. The recombination condition ( |14D is equivalent 
to 

Rij < Rpij with p i:j = — — T,J . (19) 

Therefore, choosing the same value for R, jets obtained with the cone algorithm ( |19"1) are 
slightly wider than those constructed with the hp algorithm (0). In order to demonstrate 
that our numerical routines work also with the kj< algorithm we have calculated some 
representative cross sections with this jet definition as well. 



3.2 Photon plus Jet Cross Sections 

Now we shall present our numerical results. We start with various cross sections for the 
7 + (n + l)-jet final state since we think that these cross sections, although smaller than 
the fully inclusive photon cross section, will be measured first due to reduced background 
problems. It is clear that in NLO the final state may consist of two or three jets where 
one jet is always a photon jet. The remnant jet is not counted since it is produced with 
zero transverse momentum. The three-jet sample, equivalent to 7 + (2 + l)-jets in the 
notation of the previous sections, consists of all 7 + (2 + l)-parton-level jets which do not 
fulfill the cone constraint R^j < R with Rij given in (0). The 7 + (1 + l)-jet sample 
contains events where two partons (possibly a photon) are recombined into one jet or one 
parton does not obey the cut on transverse momenta (0). In the following we shall sum 
over the two samples with n — 1 and n = 2 jets. If there are two jets in an event, we order 
them according to their transverse momenta and call the one with larger p? "jet 1" and 
the one with smaller p? accordingly "jet 2". Also, from now on we will use a simplified 
notation and denote kinematic variables of the jet containing the photon simply by Pt,-j, 
t] 7 and 7 . 

In Figs. f| and [| we show results for the px and r\ dependence of the cross sections da/dpr 
and da/dr] for the photon jet and the jet with the largest p?- In each of the four figures 
we have plotted three curves for three choices of z cut defined in (0), z cut = 0.5, 0.7 and 
0.9. Together with the cone radius R, z cnt controls the amount of photon isolation. As 
to be expected the cross section decreases with the degree of photon isolation, i.e. with 
increasing z cut . Specifically in Fig. [|a we present dcr/dpx^ as a function of px n for the 
three z cnt values and for px n > 5 GeV. All other variables, in particular r/j e t, ?7 7 and 
Prjet are integrated over the kinematically allowed ranges. We see that all three cross 
sections have a similar shape. In Fig. ||a da/dpx jeti as a function of Prjcti is shown. The 
qualitative behaviour of the cross sections for the different ,2 cu t's is similar as in Fig. |]a. 
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Figure 4: p^- (a) and rj- (b) distributions of the photon jet for z cnt = 0.9, 0.7 and 0.5 
(full, dashed and dotted curves). 
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Figure 5: pr- (a) and rj- (b) distributions of the most energetic jet for z cut = 0.9, 0.7 and 
0.5 (full, dashed and dotted curves). 
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Figure 6: Distribution with respect to the photon-jet imbalance parameter z 7 i (a) and 
the jet-jet imbalance parameter z 12 (b) for z cut = 0.9, 0.7 and 0.5 (full, dashed and dotted 
curves) . 



For the rj distributions, shown in Figs. [|b, [5]b for the photon jet and the most energetic 
jet, we have integrated over p T<1 > 5 GeV and Prjeti > 6 GeV. The shapes of the three 
curves for the different z cnt values are similar for both the cross section da/dr) y in Fig. |]b 
and the cross section da/drjj et i in Fig. |]b. We note that the 7/jeti distribution peaks at 
somewhat smaller rapidities than da/dr]^. 

Distributions with respect to pt and rj of the second jet in 7 + (2 + l)-jet events do not 
depend on the isolation cut since in this case each parton (photon) constitutes a jet by 
its own. Therefore we show corresponding figures in the next subsection when we will 
discuss the influence of the cone size R on the cross sections. 

In addition to predicting distributions in the transverse momenta and the rapidities of 
the photon and the hadron jets we also have calculated distributions for variables which 
characterize the correlation of two jets. One of these variables is 

^ = -^#^. (20) 

Note that is defined with the help of the transverse momentum of the photon jet, i.e. 
Pr i7 may include a contribution from accompanying hadronic energy. The dependence 
of the cross section on z 7 i characterizes the imbalance in transverse momentum of the 
photon and the most energetic jet. Similar variables have been used before for studies in 
the case of photon plus charm jet final states in pp collisions p4j and of two-jet production 
in ep scattering [25]. The result for da/dz^i is shown in Fig. for three photon isolation 
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cuts z cu t = 0.5, 0.7 and 0.9. The cuts on transverse momenta and the cone parameters 
are as defined before. 

For two-body processes such as the LO Compton subprocess 7*g — > jq, the final photon 
and the jet have balancing transverse momenta and the distribution is a 5-function in 
(1 — Zji). Also the fragmentation process contributes only to z 7 i = 1 since in this case 
the transverse momentum of accompanying hadronic energy is collinear with the photon, 
resulting in p?^ = PT,bare 7 + Pr,had = Prjeti- Contributions with z n \ ^ 1 are due to the 
higher-order three-body contributions. Events with z 7l < 1 typically result from config- 
urations where a single photon is opposite in transverse momentum to a jet consisting 
of two partons. Since according to the recombination prescription (|15|) the scalar sum 
of transverse momenta is ascribed to the jet, not the vectorial sum, one finds always 
Pt,7 < PTjeti and thus z 7 i < 1 in this case. Moreover, the photon is never accompanied 
by a hadronic parton in this case and events with z 7 i < 1 are consequently not affected 
by the isolation cut. On the other hand, events with z 7 i > 1 are predominantly due to 
configurations with a photon-jet consisting of a photon and a quark or gluon opposite to 
a jet consisting of a single parton. In this case, a variation of the photon isolation cut z cut 
has a strong effect on the differential cross section. 

These features are clearly visible in Fig. ^]a. The cross section increases when lowering 
z cut , i-e. when larger fragmentation contributions are included. We notice that the z 7 i 
distribution is not symmetric around z 7 i = 1. The cross section for z 7l < 1 is larger 
than for z 7l > 1, becoming more and more symmetric for less restrictive isolation cuts on 
z 7 . The residual asymmetric behaviour of this distribution for vanishing isolation cut is 
a dynamical property of the underlying cross section. 

In the region of z 7 i near unity, one of the two final state partons of three-body contribu- 
tions (not the photon) becomes soft and thus this region is sensitive to soft-gluon effects. 
In our calculation with an invariant mass cut slicing parameter y^ these soft-gluon correc- 
tions are considered as two-body contributions. Their contribution depends on the slicing 
parameter y$ . To remove this dependence, i.e. to remove the infrared sensitivity, we must 
include a sufficiently large fraction of the three-body contribution from outside z 7 i = 1. 
Therefore we averaged the z^\ distribution over sufficiently large bins and studied the sum 
of the 7 + (1 + 1)- and 7 + (2 + l)-jet cross sections da/dz^. We have chosen a bin width 
of Az 7 i = 0.2 around z 7 i = 1 and Az 7 i = 0.1 elsewhere. 

It is clear that the cross section outside the bin at z 7 % = 1 has a stronger scale dependence 
than inside this bin since only three-parton terms contribute. The cross section inside 
the bin at z 7 i = 1 is a genuine NLO prediction with expected reduced scale dependence. 
The scale dependence will be studied later for some other distributions. The 5-function 
behaviour at LO is of course in reality modified not only by NLO corrections, but also by 
non-perturbative effects originating from hadronization and a possible intrinsic transverse 
momentum of the initial parton. Our calculation does not include these latter effects. 
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Figure 7: (^-distribution of the photon jet for z cnt = 0.9, 0.7 and 0.5 (full, dashed and 
dotted curves). 

In Fig. Hb we show the cross section da/dz^ where the variable 
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PTJetl 



(21) 



with PTjeti > PT,jet2 measures the correlation between the two jets in 7 + (2 + l)-jet 
events. This cross section peaks at z 12 = as to be expected and decreases away from 
Z\2 = 0. The point z i2 = is the point with no second jet, i.e. the pure 7 + (1 + l)-jet 
region. This region is again infrared sensitive. Therefore we integrated here over the bin 
—0.05 < Z12 < 0.05. Outside this region we chose a bin size of Az±2 = 0.05. Note that the 
distribution shown in Fig. |^b includes the contribution from low-]>r partons with px below 
the cut in (|T2"|). We have also calculated the ^-distribution restricting to 7 + (2 + l)-jet 
events where both jets have px > 6 GeV. In this case, the distribution would extend 
to larger values of Z12 with a maximum at z^ ~ 0.6. The asymmetric behaviour of the 
curve visible in Fig. ^b, i.e. the tail at large Z12, is due to the contribution from these 
7 + (2 + l)-jet events. 

Another interesting variable might be the azimuthal angle 7 of the emitted photon. We 
define 7 with respect to the plane spanned by the momenta of the beam and of the 
scattered lepton. In Fig. |7| we show the dependence of the cross section on 7 , again for 
the three z cut values 0.5, 0.7 and 0.9. As is seen in this figure, the photon is emitted 
dominantly at 7 = 0. We note that the distribution becomes flatter with decreasing 
z cut . It is symmetric within the statistical accuracy of the calculation. We do not present 
here a similar plot for da/d(pj e ti, the cross section with respect of the azimuthal angle 
of the most energetic jet. It would show a distribution which peaks at 0j et i = n, since 
the dominant contribution to the cross section comes from configurations in which the 
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Figure 8: px- (a) and r\- (b) distributions of the photon jet for R = 1.0, 0.7 and 0.5 (full, 
dashed and dotted curves). 



photon and the jet with the largest pt are emitted back-to-back. For jet 2 there is no 
such correlation. The cross section da/dcj)^ is independent of 0j et 2 and, in NLO, does 
not change with z cnt - In a similar way one can discuss the cross section as a function of 
7 i = 7 — 0jcti or 7 2 = 4>-y — 0jet2- da/dtp^i is strongly peaking at </> 7 i = tt. Here it is 
again necessary to calculate da/dcf)^ with a sufficiently large bin size around 7 i = 7r in 
order to avoid any infrared sensitivity. On the other hand, the distribution da/dcf)^ is 
flat around 72 = and decreases towards 72 = n. 



3.3 Cone Size Dependence of Jet Cross Sections 

So far we presented results only for the cone jet algorithm with the cone radius fixed to 
R = 1 for both the photon jet and purely hadronic jets. Sometimes it is advantageous 
to use smaller cone radii to suppress background processes. On the other hand the 
dependence of the cross sections on the cone radius is a genuine NLO effect since LO 
cross sections do not depend on the jet definition. To present an overview of the cone 
radius dependence we have recalculated some of the cross sections shown so far for R = 1 
for two smaller radii R = 0.5 and 0.7. For the photon isolation parameter we fix now 
z cut = 0.9. 

In Figs. |8]a, b we show the results for da/dpT^ and da/drj^ with R = 0.5, 0.7, and 1.0. All 
other cuts are chosen as before. The case with R = 1 and z cut = 0.9 was shown in Figs. 
f|a, b already. The distributions for the jet with largest pr are exhibited in Figs. |9]a and 
b. da/dpr^, da/drj^, da/dprjeti and da/dr]- ]et i show very little dependence on the cone 
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Figure 9: px- (a) and r\- (b) distributions of the most energetic jet for R = 1.0, 0.7 and 
0.5 (full, dashed and dotted curves). 



size R. For a less restrictive isolation cut of z^ > 0.5 these distributions would decrease 
with decreasing cone radius R as it is known from other jet calculations. 

In our previous work |TT| we studied the equivalent cross sections also for the two event 
classes with 7 + (1 + l)-jets and 7 + (2 + l)-jets separately. It turned out that the contri- 
bution for 7 + (1 + l)-jets is the dominant one. This is expected, since the contribution 
with 7 + (2 + l)-jets is an 0(as) effect. Here we present now also results for the pr and rj 
distributions of jet 2 in this latter process. The results are shown in Fig. [l0| a (da / 'dpr, j e t2) 
and Fig. |10| b (da j drj^) ■ For £>T,jet2 > 6 GeV the cross section is very much reduced as 
compared to the cross section in Fig. |9|a. In fact, the dominating event configuration is 
with a photon and one jet balancing each other in p^', a third jet with comparable px is 
found in only a small portion of the events. The rapidity distribution da/dr]j Ct2 plotted in 
Fig. |T0"| p has a larger tail extending to larger rapidities as compared to the cross sections 
da/dr] 1 and da/dr]^. The second energetic jet originates dominantly from the incom- 
ing quark and therefore is in many cases closer to the proton remnant, i.e. at positive 
rapidities, than the harder jet. 

The cross sections da/dpTj e t2 and da/drj )et 2 increase with decreasing R. This is the be- 
haviour expected for cross sections which are of leading order in a s . 

In our previous work |L2] we studied da jdz^ as a function of the fraction z 1 of the mo- 
mentum of the photon inside the photon jet. This cross section is expected to contain 
information on the photon fragmentation functions. The results presented in [12 were 
obtained for the case R = 1 only. To see how results change with R, we show in Fig. |TT 
da j z 1 for the particular photon fragmentation function of Bourhis et al. [19] for R = 0.5, 
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Figure 10: (a) and r\ (b) distributions of the second, less energetic jet for R = 1.0, 0.7 
and 0.5 (full, dashed and dotted curves). 



0.7 and 1.0. As is seen, da/dzj decreases with decreasing R since this cross section is a 
superposition of leading and next-to-leading order contributions. 
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Figure 11: distributions for R = 1.0, 0.7 and 0.5 (full, dashed and dotted curves). 
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Figure 12: p?- (a) and r\- (b) distributions of the photon jet for R = 1.0, 0.7 and 0.5 (full, 
dashed and dotted curves) in the kx algorithm. 



The R dependence of the other cross sections considered above follows the same pattern. 
In cases where we have a superposition of LO and NLO contributions we encounter a 
decreasing cross section with decreasing R. In regions where the cross section receives 
contribution from 0(a s ) only with no additional NLO corrections included, the cross 
section increases with decreasing R. Thus for example da/dz 7 i (cf. Fig. ||a) is decreasing 
with R in the bin at = 1 and increasing with R outside this bin, i.e. for z 7 i < 0.9 and 
> 1.1. The R dependence of dajdzyi (cf. Fig. ||b) is similar: inside the bin around 
zyi = the cross section decreases with decreasing R and outside zu = it increases with 
decreasing R. 



3.4 Jet Algorithm Dependence 

Measurements of cross sections for the production of a photon plus jets, as for example 
in 7p collisions ||, have been performed with the help of the cone algorithm used to 
define jets and the photon isolation. However, the kx algorithm has definite advantages, 
in particular in the experimental analysis [F22J. Therefore we present a few distributions 
based on the k? algorithm as well. This algorithm is used for the recombination of two 
partons into a jet as explained in sect. 3.1 as well as for the definition of the photon jet. 
In Figs. |l^a, p~3| a, we show da/dpr^ and da / 'rfprjeti calculated with the kx algorithm for 



R = 0.5, 0.7 and 1. These two groups of curves have to be compared to the results with 
the cone algorithm in Figs. [8]a, §a. The qualitative behaviour of the two cross sections 
is similar; but there are quantitative differences. We note that the dependence on the 
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Figure 13: px- (a) and r\- (b) distributions of the most energetic jet for R = 1.0, 0.7 and 
0.5 (full, dashed and dotted curves) in the kx algorithm. 



parameter R, which controls the size of the jets, is now even more reduced as compared 
to the corresponding cross sections with the cone algorithm. The cross sections da/drj^ 
and da/dr]j et i for the kx algorithm are displayed in Figs. 12b, 13 b. They can be compared 
with the corresponding cross section for the cone algorithm in Figs. |8|b, |S|b. Similar to 
the px distributions the qualitative behaviour did not change. Again the R dependence 
seems to be reduced for the kx algorithm. This is even more the case for the cross 
sections da/dpx,] e t2 and da/dr]j et 2, shown in Figs. |T4|a, b which should be compared to the 

b. It is clear that the cross sections in Figs. 14a, b 



corresponding results in Figs. |10|e 
increase with decreasing R in the same way as the cross sections for the cone algorithm 
in Figs. |l0la, b. 

A direct comparison of cross sections calculated with either the cone algorithm or the kx 
algorithm is shown in Fig. [15] (dashed and dotted curves). Here we have chosen R = 1 and 
-2 C ut = 0.9. In Figs. |l5|a, b the cross sections da/dpx n and da/drj^ are plotted, respectively. 
As we can see, these cross sections hardly change when the cone algorithm is replaced by 
the kx algorithm. Only where da/drj^ is maximal, i.e. near 7]^ = —1.5, the cross section 
with the kx algorithm is approximately 5 % larger than with the cone algorithm. 
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Figure 14: p T - (a) and redistributions (b) of the second jet for R = 1.0, 0.7 and 0.5 (full, 
dashed and dotted curves) in the kx algorithm. 




Figure 15: Comparison of the p?- (a) and ^-distributions (b) of the photon jet in the cone 
algorithm with z 7 > 0.9 and R = 1.0 (dotted curves), in the kx algorithm with R = 1.0 
(dashed curves) and without jet analysis (inclusive cross section, full curves). 



22 



It is expected that the requirement to observe additional jets reduces the cross section 
for the production of a high-p^ photon. To study this reduction we have calculated also 
da/dpT,-y and da/drj^ for the inclusive case, i.e. without additional jets required in the 
final state and the same photon isolation constraint in terms of an isolation cone around 
the photon and the cut z cut = 0.9. The results for these inclusive cross sections are shown 
in Figs. |T5|a, b as full curves. The total cross section without jet algorithm applied is 
increased by about 35%. This is reflected in the inclusive cross section da/drj^ which is 
always larger than the corresponding distribution for final states which are required to 
contain at least one jet, over the full range of r] 1 . In da/dpx^ only the first two pr-bins 
contain an appreciably larger cross section due to the removal of the jet requirement. 
This can be traced back to the different prescriptions used to remove events with small 
transverse momentum: for the analyses based on the cone or the jet algorithms, we 



used the cuts Q12Q which are applied to individual jets, whereas in the case without 
jet algorithm we applied the cut ( JT~5| ) to the sum of transverse momenta of all hadronic 
particles in the final state. As a consequence, an event with two low-pr partons, each with 
Pt < 6 GeV will be rejected in the first case if these two partons are not recombined into 
one jet, while the event will be accepted in the second case if the sum of the transverse 
momenta of the two partons is larger than the cut of 6 GeV. This clearly affects only the 
bins with lowest pt- 



3.5 Scale Dependence of Jet Cross Sections 

All results presented so far have been obtained for a renormalization and factorization scale 
fj, fixed at /i = prjeti (/-t = Et of (pT3|) for inclusive cross sections). In LO cross sections, the 
scale dependence is exclusively due to variations of /i in the parton distribution functions. 
At NLO we expect that additional terms containing an explicit /x-dependence will reduce 
the scale dependence. Instead of studying the scale dependence of all the differential cross 
sections discussed so far separately, we have investigated the scale dependence of some 
components of the total cross section, i.e. integrated over the phase space allowed by the 
transverse momentum cuts (|l^). We define the scale in the form jj 2 = f 2 p\- ]Gtl and vary 
/ between / = 1/4 and 4. 

In Fig. |16| a we have plotted the / dependence of the 7 + (1 + l)-jet cross section in 
LO and NLO (denoted 0(a s )) and of the 7 + (2 + l)-jet cross section. The LO cross 
section (denoted "Born" in Fig. ^) increases with / by approximately 10 % in the range 

0. 25 < / < 4. The cross section including corrections of 0(a s ) is almost independent of /, 

1. e. no scale dependence inside the considered range of scales is visible. Here the decrease of 
the cross section due to the decrease of a s with increasing / is compensated by the increase 
originating from the scale dependence of the proton PDF's. The cross section for the 7 + 
(2 + l)-jet final state, being of order 0(a s ), decreases with increasing / by approximately 
25 %. This is a combined effect of the dependence of a s and of the / dependence of the 
parton distribution functions. The / dependence of separate components ("real", "sing" 
and "frag") of the 7 + (1 + l)-jet cross section including 0(a s ) terms is plotted in Fig. 
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Figure 16: Scale dependence of the total cross section (a) and of separate contributions 
to the 7 + (1 + l)-jet cross section (b). 



Top and again compared with the LO cross section, "real" stands for the tree graph 
contributions in 0(a s ), calculated with the slicing parameter as described in section 2.4. 
Since it is a tree-graph term it decreases with increasing / due to the decrease of a s . The 
singular term "sing" , which includes virtual corrections and singular contributions below 
the slicing cut, is negative and decreases in absolute value by the same amount as "real". 
The fragmentation contribution "frag" which includes both terms of the right-hand side 
of (0) is also negative and almost independent of / as expected, since the factorization 
scale dependence (fip in (^)) cancels in first approximation. 



3.6 Dependence on Low-p^ Cuts 

The choice of two different cuts for the transverse momentum of the photon and the jet 
is needed to avoid the otherwise present infrared sensitivity of the NLO predictions. This 
sensitivity is known from similar calculations of dijet cross sections in ep collisions [^5[ 
and must be avoided. The same problem was encountered in the calculation of inclusive 
two-jet cross sections in 7p collisions ||26|| , for the production of a prompt photon plus a 



charm quark in pp collisions EM and much earlier in NLO calculations of the inclusive 



cross section for photon-hadron and for two-photon production f27 



Above we have chosen the difference between the two pr,min cuts for the photon and the 
jet, A = p™j — Pt™, equal to 1 GeV (see (|T2"|)). A should not be too small since then 
we would encounter the infrared sensitive region where the prediction of the cross section 
becomes unreliable. In order to obtain some information about possible choices for A we 
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Figure 17: Dependence of the total cross section on lower transverse momentum cutoffs: 
as a function of p™j with fixed p™™ = 5 GeV (a) and as a function of p™™ with fixed 
= 5 GeV (b). 



have studied the 7 + (n + l)-jet cross section da/dp TJ integrated over p T j > p™ m as a 
function of The transverse momentum of the photon was always integrated over 



the range > Pt™ = 5 GeV. The results for <j(ptj) are plotted in Fig. [T7|a. Starting 
at p™j = 6 GeV this cross section increases with decreasing p™j with almost constant 
slope. At about p™™ = 5.5 GeV the slope decreases and approaches zero and even changes 
sign so that <r(p™j) develops a maximum below p™j = 5.5 GeV. This change of slope is 
due to the infrared sensitivity in the point p™j = p™™. To avoid this region one must 
choose A^O. From the plot we observe that A > 0.5 GeV would be sufficient. Thus, in 
principle, we could have used a smaller value for this difference than was chosen in (|T2|). 

The cross section in the vicinity of p™j = p™™ cannot be predicted reliably. It depends 
on the technical cut as soon as one approaches the limit A = 0. In fact, the infrared 
sensitive region is very much influenced by non-perturbative effects which are not included 
in our calculation. In any case it would be interesting to measure a(p^ m ) in order to 
investigate this non-perturbative region. 



In Fig. |17j a we present cr(p™j) also for p™ l j below 5 GeV, i.e. for A < 0. If p™ l j is 



increased starting from p™ m = 4 GeV, the slope of cy{p™ l j) changes at about p™ j = 4.5 
GeV and crip^j) becomes smaller towards p™j = 5 GeV compared to the behaviour with 
constant slope. This stronger decrease of cr(p™j) above p™™ = 4.5 has its origin again in 
the infrared sensitivity of this region. For p™ m — > 5 GeV, and A < the cross section 
approaches the same value as we have obtained for A — >• in the region A > 0. This 
means, o"(p™j) is not singular at p™™ = 5 GeV, nor has it a discontinuity. 
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The observed behaviour is of course not only visible at the specific value of p™ m = 5 GeV 
which we have chosen. A similar dip is seen for larger and smaller values as well. It 
becomes more pronounced for smaller values and gets washed out for larger p™ m . 

For completeness we show in Fig. |TT| b also the cross section with the roles of p^j and 
interchanged, i.e. we fix p™™ = 5 GeV and study the dependence on the cut for the 
photon transverse momentum The behaviour of the cross section in Fig. |T7| b is very 
similar to the first case and shows the same infrared sensitive region. Inside the region 
—0.5 GeV < A < 0.5 GeV the cross sections agree inside numerical errors. Only for 
larger | A| the dependences on the minimal transverse momenta of the jet and the photon 
become different. The cross sections shown in Fig. [17| have been calculated with the cone 
algorithm using R = 1 and z cnt = 0.9. 



4 Summary and Concluding Remarks 

We have presented results of a next-to-leading order calculation of isolated photon pro- 
duction in large-Q 2 ep scattering. Contributions from quark-to-photon fragmentation are 
explicitly taken into account. We have discussed numerical results for 7 + (1 + l)-jet and 
7 + (2 + l)-jet cross sections as functions of transverse momenta, rapidity and other ob- 
servables derived from photon and/or jet kinematic variables. Infrared sensitive regions, 
as for example the region of equal photon and jet pr, are studied in detail. 

We investigated several of these cross sections with respect to their photon isolation, jet 
cone size, scale and jet algorithm dependences. It was found that these dependences are 
rather weak. In particular, the scale dependence of the integrated cross section is very 
small, giving quite some confidence in the reliability of our predictions. Also the results 
depended very little on the choice of modern parton distribution functions for the proton. 

We expect that the measurement of photon plus jet cross sections at HERA will contribute 
to testing perturbative QCD in the process 7*p — > 7A, in an area which has not been 
studied yet. The calculation covers the range of large Q 2 where the results do not depend 
on parton distribution functions of the virtual photon. At even larger Q 2 additional 
contributions from Z boson exchange become important, which have been neglected in 
the present calculation. It is no major problem to incorporate these missing parts. They 
are non-negligible and must be considered when experimental data become available at 
very large Q 2 . 
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